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Na 4 lr 30 g is a unique case of a hyperkagome 3D corner sharing triangular lattice which can be 
decorated with quantum spins. It has spurred a lot of theoretical interest as a spin liquid candidate. 

We present a comprehensive set of NMR data taken on both the ^®Na and sites. We show 
that disordered magnetic freezing of all Ir sites sets in below T/ ~ 7 K, well below J = 300 K, 
with a drastic slowing down of fluctuations to a static state revealed by our Ti measurements. 
Above typically 2 T/, physical properties are relevant to the spin liquid state induced by this exotic 
geometry. While the shift data shows that the susceptibility levels off below 80 K, l/Ti has little 
variation from 300 K to 2 T/. We discuss the implication of our results in the context of published 
experimental and theoretical work. 


Geometric frustration has long been known as an essen¬ 
tial ingredient to stabilize a quantum spin liquid (QSL) 
state in more than one dimension (1D)[I1[2]. Since ex¬ 
perimental realizations are rare, Herbertsmithite 0 trig¬ 
gered a lot of excitement with its discovery in 2005. It is 
the most well known example of a two dimensional (2D) 
QSL where -S' = ^ Cu^+ spins form a perfect kagome lat¬ 
tice of corner-sharing triangles with dominant Heisenberg 
interactions, no sign of ordering [HIS] and a continuum of 
excitations interpreted as a fractionalization of excita¬ 
tions into S = ^ spinons [5], similar to ID S' = ^ Heisen¬ 
berg chains. Recently other routes to 2D QSL physics 
have been explored such as frustration induced by com¬ 
peting interactions with Kapellasite mi and, outside 
of cuprates, the kagome vanadate DQVOF with S = | 
[5], Mo®+ double perovskite (TD], and organic com¬ 
pounds which are driven by proximity to a Mott transi¬ 
tion nuiin. 

Currently, Na4lr308 is one of the most compelling 
frustrated QSL candidate in three-dimensions (3D) [TSj 
where Ir'^"*' ions with effective Jeff = \ form a corner¬ 
sharing lattice of triangles named ’’hyperkagome” El- 
Strong spin-orbit coupling, SOC, has been identified as 
an important ingredient in the Hamiltonian. Several pos¬ 
sible theoretical scenarios have been proposed thus far, 
including a classical long-ranged order-by-disorder 120° 
coplanar ground state m which in the quantum limit 
melts into a gapped QSL [1^ and a gapless QSL [TTIIT^ . 
More generally, iridates appear as an ideal playground 
for the study of novel physics governed by strong SOC 
in competition with Coulomb repulsion, crystal field ef¬ 
fects, and inter-site hopping. This has led theorists to 
promote, for example, the Heisenberg-Kitaev model m 
and the spin-orbit Mott insulator [20]. Na4lr308 also ap¬ 
pears to be close to an insulator-metal transition [2TH25] 
and weak Mottness has been proposed as a possible sce¬ 
nario for a spin-liquid ground state muz], like in the 
case of organic triangular spin-liquids El- 

The macroscopic susceptibility, X) is typical of Jeff = | 
moments with antiferromagnetic interactions J ^ 300 K 


EIIIH]- Heat capacity, C, shows no sign of a bulk transi¬ 
tion and has two remarkable features: at 24 K ~ J/10 a 
broad maximum and at low temperatures C = qT + /3T" 
where 2 < n < 3 [a El]. Both X and C suggest the pres¬ 
ence of a gapless ground state [niiiH]- Magnetocaloric 
measurements suggest a quantum critical behavior in 
zero-field m- Although no signs of a bulk transition 
exist, there is a small splitting of the field-cooled (FC) 
and zero field-cooled (ZFC) magnetization at T ~ 6 K 
P5] . Initially, this splitting was proposed to be associ¬ 
ated with a ~ 1% defect or impurity term [THIH]: but 
recent /rSR measurements EH suggest quasi-static short 
range spin correlations appear below T = 6 K. 

Probing the role of these defects on the physics of this 
compound, disentangling their contribution to the sus¬ 
ceptibility from that of intrinsic origin and revealing the 
low energy spin dynamics are the central unexplored fo¬ 
cus of the ^^Na and NMR studies presented in this 
Letter, down to 1.2 K . After explaining the site assign¬ 
ment and respective advantages of both probes, we are 
able to study different aspects of the underlying physics. 
The shift points to a plateau in the local suscep¬ 
tibility at low temperatures which may be viewed as a 
partial confirmation of the predictions of HUES]- The 
evolution of the broadening of and ^^Na lines reveals 
a bulk transition to a frozen state at T ~ 7 K with a 
magnetic moment, fiir 0.27fj,B- Relaxation data, 
on ^^Na line confirms a phase transition at T = 7.5 K 
while at high temperature, it has very little T-variation 
which serves as an additional signature of the spin liquid 
regime in this exotic geometry. 

The synthesis and the quality assessment of our sam¬ 
ple are detailed in [3^. NMR measurements on ^^Na, 
nuclear spin J = | and gyromagnetic ratio ^^7/27r 
= 11.262 MHz/T, and / = § and 1 ^ 7/271 = 

5.772 MHz/T, were performed in a pure sample and, for 
comparison, in a depleted sample (Na 4 _a;Ir 308 ), in fixed 
and variable field magnets in a range of 4 T to 11 T us¬ 
ing home-built probes. Fourier-transformed spectra were 
obtained using a standard Hahn echo pulse sequence. 
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FIG. 1. Spectra taken at 300 K with a frequency of 68MHz 
for ^®Na in Na4lr308 (i) and Na4_a;Ir308 (ii) and with a fre¬ 
quency of 43.3MHz for in Na4lr308 (v). Solid lines rep¬ 
resent raw data; Gaussian curves indicate the location of the 
different sites: peak A as blue, peak B as pale blue, peak C 
as light gray, peak D as dark grey, peak E as black. The local 
environment of (hi) 0(1) and (iv) 0(2) are shown. 


The frequencies are given in comparison to a ^^Na and 
reference, NaCl solution, used as the zero frequency. 
Spin-lattice relaxation rates were obtained with a satu¬ 
ration recovery pulse sequence. 

The crystal structure has three ^^Na sites. Na(l) 
combined with three Ir form the tetrahedron of the py- 
rochlore lattice; Na(2) and Na(3) also form a network of 
corner-shared tetrahedra, each having a 75% occupancy 
M- In the unit cell, the ratio of Na( 1/2/3) is calcu¬ 
lated as (4:3:9). The spectrum of the central transition 
of Na4lr308 is shown in Fig. [^i). While the three ex¬ 
pected sites could not be resolved, the spectrum has two 
distinct peaks, labeled A and B, with weights of 75(6)% 
and 25(6)% when corrected for the relaxation rates, re¬ 
spectively. Since Na(3) is 56% of the total Na, part of 
peak A must be associated with Na(3). A further anal¬ 
ysis of the first order satellite lines [SD] supports this 
conclusion. Comparing the spectrum for Na 4 _a;Ir 3 08 
to Na4lr308 (Fig. [^ii) to (i), respectively), peak B has 
lost a considerable amount of intensity. Since the occu¬ 
pancy of Na(2) and Na(3) decreases with increasing x in 
Na4_2,Ir308 [Mj, peak B is most likely associated with 
the Na(2) crystallographic site. The last site, Na(l), was 
impossible to resolve but is expected to account for the 
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FIG. 2. The shift of the three sites (markers, left axis) 
and the bulk susceptibility (solid line, right axis) as a function 
of temperature. Inset: (i) spectra of peaks G and D 
taken with a frequency of 68.077 MHz from 300 K to 100 K 
in 40 K intervals, (ii) Isolated spectrum of peak E obtained 
by contrast with a frequency of 37.518 MHz from 100 K to 
20 K in 20 K intervals and 15 K. The blue markers show the 
bulk susceptibility. 


rest of the intensity of peak A. 

has two crystallographic sites, 0(1/2), with ratios 
1:3 whose local environments are shown in Fig. [^iii) and 
(iv), respectively. The spectrum is shown in Fig. with 
three peaks labeled C, D, E. Since Na(2/3) have a 75% 
partial occupancy [14], the 0(1) site has two dominant 
local environments which, assuming a random binomial 
occupation, should each account for 10.5% of the total 
spectral weight: one fully occupied with three-fold rota¬ 
tional symmetry and one with two of three Na sites oc¬ 
cupied with no rotational symmetry. An analysis of the 
first order satellites [30] shows that peak C is the only 
site with high symmetry and therefore corresponds with 
0(1). Since peak D has a similar spectral weight to peak 
C (17(5)% and 15(3)% respectively), it also corresponds 
with 0(1). The remaining peak E which has the bulk of 
the spectral weight (67(8)%) must correspond with 0(2). 

Most of our NMR static data reported in this Letter 
have been taken on O sites which are better coupled to 
Ir moments. We used a combination of the data from 
the 0(C,D,E) lines in order to analyze the shift and the 
linewidth, our selection dictated by resolution and cou¬ 
pling considerations. The temperature evolution of the 
spectra are an example (Fig. insets (i) and (ii)). The 
0 (1,0) peak is resolved the best and can be followed 
down to 20 K. The 0(2,E) peak can be isolated by a 
contrast experiment using fast repetition rates. Its shift 
can only be reliably extracted down to 40 K due to sub¬ 
stantial broadening. 

In the paramagnetic regime, the spectral shift is lin¬ 
early related to the local susceptibility, xioc, by K{T) = 
Kq -b AhfXioc(T’) where Kq is a temperature independent 
orbital contribution and Ahf is the hyperfine coupling 
constant [3D]. The spectral shift tracks the intrinsic sus- 
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FIG. 3. FWHM of the Gaussian broadening with respect to 
15 K in kOe for ^^Na at 4 T and 7 T (left axis) and at 
7 T (right axis) as a function of temperature below 16K. The 
fj,SK data from [55], normalized to the NMR IK AH, as pink 
triangles. Inset: ^^Na spectrum at 45.046 MHz for 4 T (solid 
line) and 78.937 MHz for 7 T (empty circles) at T = 1.3 K. 
The spectrum centered at 4 T has been shifted horizontally 
by 3.005 T and renormalized for comparison. Background 
Gu is from the experimental setup. The blue line indicates 
the expected shape of the spectrum if long ranged order were 
present with the same value of the moment. 

ceptibility while the magnetic broadening of the line re¬ 
flects the distribution of local susceptibilities, in general 
due to a small amount of defects. Thus, spectral mea¬ 
surements can uniquely disentangle these two contribu¬ 
tions to the macroscopic susceptibility. We find a good 
agreement between K for 0(1,C), 0(1,D) and 0(2,E) 
peaks and Xmacro between 300 and 100 K and establish 
that the intrinsic susceptibility levels off below 80 K. The 
lower T-limit of 20 K of our data conservatively indicates 
the absence of a gap larger than J/15. 

Looking closer at Fig. [^ii), the linewidth, AH, of the 
0(2,E) peak gradually increases down to 15 K as ex¬ 
pected from the Curie tail of Xmacro- Below 10 K, AH 
increases significantly more, as seen in Fig.[^ This broad¬ 
ening is best tracked by the ^^Na(A) line which remains 
well resolved, has a small high-T width and can be stud¬ 
ied at various fields due to the Na NMR sensitivity m- 
Below 7.5 K where the ZFC-FC susceptibility measure¬ 
ments split, ^^AH is the same for 4 T and 7 T which 
is a landmark of static freezing. This frozen phase is 
fully settled by 2K below which the linewidth is con¬ 
stant, consistent with /iSR measurements [55]. If this 
broadening were due to the onset of an ordered phase, 
the well-defined internal fields would cause the powder 
spectrum to be rectangular. As shown in the inset of 
Fig. § this is not the case and the identical spectra at 
1.3 K can be approximately fitted by a Gaussian indi¬ 
cating a disordered freezing of Ir moments. The absence 


of any residual narrow spectral contribution which would 
be typical of a paramagnetic phase indicates a freezing 
of all moments. 

In the frozen phase, the FWHM of the Gaussian broad¬ 
ening with respect to 15 K tracks the order parame¬ 
ter as in the muon data [29] . Since the hyperfine con¬ 
stant 26(5) kOe/^_B is most accurately determined for 
the 0(2,E) peak from the K — x plot in the T > 100 K 
paramagnetic regime (see Fig. [^, the magnetic moment 
of Ir^'*' can be best extracted from at 1.3 K . In 

the framework of a Gaussian lineshape, /iir = - F — 
where z is the number of Ir coupled to the O. This yields 
0.27(4 )^b. This estimate is comparable to the observed 
/Tir from neutron measurements in other iridates such as 
Sr 2 lr 04 [SlEg and NasIrOg [SS]. 

We now turn to the spin dynamics as revealed by our 
relaxation data. We have measured on the ^^Na(B) 
peak which is well-defined until below 15 K. It has suffi¬ 
cient coupling to the Ir site with better separation than 
the O sites which allows us to isolate the relaxation com¬ 
ponents. (23) For the nuclear spin of ^^Na, /=§, the 
magnetization recovery is described by only two compo¬ 
nents [33] 

M(t) = Mo[l - + (1 - a)e-(*/^il'')] (1) 

where Mq is the equilibrium magnetization, a measures 
the weight of the short and long components of Ti and 
is typically /3 1 is introduced here to account for 

the distribution of relaxation rates [35] which is related 
to the degree of disorder in the system. 

decreases very smoothly until T ~ 24 ^ with a 
phenomenological fit = A + BT°‘ (Fig. where 
the constant term A = 0.010(2) ms“^ dominates. The 
300 K value is consistent with the exchange narrowing 
limit of Heisenberg antiferromagnets given by Moriya’s 
theory [36], 1/Ti ^ A^/J = 0.018(9) ms“^. Given the 
only slight decrease of 1/Ti, 1/TiT increases when T de¬ 
creases which could indicate a moderate strengthening of 
correlations until reaching a maximum at the phase tran¬ 
sition at T ^ 7.5 K. Below 7.5 K, drops by two orders 
of magnitude which is typical of the critical slowing down 
of spin fluctuations when entering a frozen phase. The 
relaxation rates become distributed below 30 K where 
j3 begins to decrease monotonically (Fig. Qi))- For a 
gapped state, ^ exp(—A/T) at low temperatures. 
An exponential (A ~ 12(2) K) and a phenomenological 
power law (power ~ 3.8(3)) fit were made for T <7 K. 
There is no substantial difference between these two fits. 
Both fits require a vertical offset of ^ 0.001 which might 
originate from a minute amount of paramagnetic uncou¬ 
pled moments. We note that the absence of gapped be¬ 
havior in the specific heat might favor the power law for 
1/Ti. Quite surprising is the discrepancy between our 
l/Ti.NMR spin lattice relaxation rate and the ^SR relax¬ 
ation rate I/Tl^^sr [IS]. l/Ff^^iSR keeps constant which 
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FIG. 4. ^Sy-i axis, black circles) and C/T from [TT] 
(right axis, blue squares) as a function of temperature. The 
purple lines are fits. Inset: (i) The stretch exponent, j3, as 
described in the text as a function of temperature, (ii) The 
relaxation curve (solid circles) at 1.3 K fit with a stretched 
exponent (blue line, /3 = 0.4) and exponential fit (grey line, 

/3 = 1). 

was argued to be the signature of configurationally de¬ 
generate phases with fluctuating order [29] . Recent LDA 
calculations on model systems I371I3HI show the /r’*', as a 
positive charge, might induce local deformations in its en¬ 
vironment which could affect the dynamics of Ir^+ close 
to the muon probe. 

From the experimental results to be discussed below, 
two T-regions emerge from our study: one typical of the 
ideal hyperkagome, Na4lr308, at high T and one whose 
physics deviates from the ideal Hamiltonian(s) proposed, 
leading to a fully frozen phase below ~ 7 K. 

Based on our NMR data we draw the conclusion that 
a hulk and static disordered spin freezing occurs at 
with inhomogeneous dynamics in the ground state. Bulk 
spin freezing might appear to conflict with the signature 
of a ZFC-FC splitting in Xmacro which was attributed 
to ’’defect” spins. This could have been interpreted as a 
marginal spin glass transition implying a weak coupling 
of such defects, e.g. mediated through a spin liquid back¬ 
ground but these defect spins, which dominate the low-T 
susceptibility, might act simply as a fingerprint of the 
bulk physics. 

In view of the 75% partial occupancy of Na [T3], we 
propose that some local charge disorder occurs. This cer¬ 
tainly affects the Ir^+ environments and leads to a dis¬ 
tribution of magnetic interactions on the hyperkagome 
lattice, especially if interactions are dominantly due to 
direct exchange as argued in many papers [33] ■ Devia¬ 
tions to equilateral interaction triangles have always led 
to transitions at T ^ J such as observed in Volborthite 


[TO] and Vesignieite [nma kagome-based compounds. 
Digging out why a distribution of interactions could lead 
to such a weak or no signature in x and C/T respec¬ 
tively should be addressed theoretically which may help 
discriminate between models proposed for this iridate. 
We note that in all investigated models, the and pos¬ 
sibly T-linear ultimate behavior of the specific heat has 
played a central role. This should be now toned down 
in view of the freezing evident in our data which inval¬ 
idates the use of T < Tf thermodynamic quantities as 
characteristic of the spin liquid behavior. 

Above the freezing temperature, a broad range of tem¬ 
peratures probes the effects of frustration on the spin liq¬ 
uid state since J/2T/ ^ 20. First, in the Tf-2 Tf range, 
the increase of 1/Ti may be interpreted as a slowing down 
of magnetic fluctuations. Whether this extended range of 
critical fluctuations could be explained by disorder or it 
rather appears as a signature of a crossover region coin¬ 
ciding with the maximum of the specific heat at T ~ 3 Ty 
remains speculative. Above ST/, three landmarks of the 
spin liquid regime now clearly appear, the pseudo-gap 
like behavior of Tfl^, the leveling off of y at 80 K con¬ 
firmed by our shift measurements, which both add up 
to the broad maximum in C/T at 24 K. Various models 
have been explored which we discuss with respect to our 
results: 

(i) A fermionic approach naturally leading to a spinon 
Fermi surface. The maximum of C/T, far too high for 
a spin glass freezing [33], could be the landmark of a 
crossover from a U{1) spin liquid to a Z 2 one with a 
paired spinon state and line nodes in the gap below 20 K 
[18] . The mixing with triplet states induced by SOC or 
Dzyaloshinkii-Moriya interactions could explain why the 
susceptibility keeps its Pauli-like behavior. Yet, our 1/Ti 
data, if intrinsic, contradicts the existence of a gap as it 
does not decrease, below the maximum of the specific 
heat. Furthermore, for a Dirac U(l) spin liquid, one ex¬ 
pects 1/Ti ~ T’* where rj is related to the shape of the 
correlation function and remains unknown [44] . This is 
not what is observed here since Ti is dominated by a 
constant term for T > 20 K. 

(ii) The transition to metallicity, observed either under 
pressure [33] or in depleted Na samples [HHH], might 
indicate the proximity to a quantum critical point. In 
this context, modeling has focused on the metallic rather 
than insulating side [MlEl]; here, the T—dependence of 
Tfl^ is a crucial test. 

(hi) A 72-sites valence bond crystal has also been pro¬ 
posed as a trial ground state [46] but, with a gap of the 
order of J, is not relevant. However, a valence bond 
glass state with a transition at ~ J/IO as argued for the 
kagome lattice is worth further exploration [37] . The low- 
T behavior of C{T) might originate from the free energy 
landscape typical of disordered systems such as conven¬ 
tional spin glasses [33] . 

In conclusion, our Ti data opens up the space of phys- 
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ical quantities to be discussed in future advanced model¬ 
ings altogether with the constant susceptibility and the 
maximum of the specific heat which have been a central 
issue in the to-date discussions. Whether a slight disor¬ 
der in the interactions might bear an explanation for the 
disordered static freezing with no signature in the spe¬ 
cific heat is also a major avenue for future work on the 
hyperkagome lattice. 
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